APS/123-QED 



The Role of Nonlinear Friction in the Dewetting of Thin Film Polymers 
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The study of the dewetting of very thin polymer films has recently revealed many unexpected 
features (e.g. unusual rim morphologies and front velocities) which have been the focus of several 
theoretical models. Surprisingly, the most striking feature of all, that is a decrease of the rim width 
with time, have not yet been explained. In the present letter, we show how the combined effects of 
a non-linear friction between the film and the substrate, and the presence of residual stresses within 
the film, result in the presence of a maximum in the time evolution of the rim width. In addition, 
we show how the introduction of a non-linear friction can also simply explain the rapid decrease of 
the dewetting velocity with time observed experimentally. 

PACS numbers: 68.60.-p, 68.15.+e, 61.41. +e, 83.50.-v 
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Since the use of nano-devices is nowadays getting more 
and more important, the stability of thin polymer films 
has become a major scientific and technologic research 
area 0]. In particular, the dewetting of thin polystyrene 
(PS) films has been extensively investigated over the last 
decade 0. A great interest has been payed to the un- 
usual rim morphology 0, 0, first studied by Reiter [(| . 
A detailed experimental study of the dewetting dynamic 
has been published two years ago 0, showing a sharp 
decrease of the dewetting velocity with time (V ~ t^ 1 ). 

These observed features have been recently explained 
theoretically in [j| . The principal ingredients allowing to 
understand the observed asymmetric shape of the rim are 
the weak friction of the film onto the substrate (grafted 
or absorbed liquid polydimethylsiloxane (PDMS) mono- 
layer), allowing a plug- flow description, and the ineffec- 
tiveness of the surface tension of the liquid at this stage 
of the dewetting process. The friction dumps the sliding 
velocity within the film over a distance A = y/ T)ho / Co 
(where r\ is the viscosity of the liquid, assumed to be 
Newtonian, hp is the initial thickness of the film, and 
Co is the friction coefficient between the film and the 
substrate) 0, leading to an asymmetric rim of width 
~ A, whereas the height of the rim H increases pro- 
gressively. Omitting inertia, the dewetting dynamic can 
be simply solved by balancing the capillary power by 
the dissipation due to friction: \S\V ~ (oAV 2 (where 
S is the spreading parameter |lOj, assumed to be neg- 
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FIG. 1: Film geometry : h(r,t) is the profile of the film, ho 
is the initial height of the film, H (t) is the height of the front, 
L(t) is the dewetted distance, W(t) is the width of the rim, 
and v(x, t) is the velocity of the film. 



ative). For a Newtonian liquid, this balance gives a 
constant dewetting velocity V — \S\/( A ,9j. The de- 
crease with time of the velocity observed experimentally 
is explained by the viscoelastic properties of the poly- 
mer film, i.e. by the fact that the liquid behaves like a 
low viscosity (770) Newtonian liquid at very short times 
(i < To = rjo/G), like an elastic solid (of elastic modulus 
G) in an intermediate time regime (tq < t < t\ = rji/G), 
and like a very viscous Newtonian liquid (with a viscos- 
ity 771 3> 770 , where 771 is related to the molecular weight 
of the polymer) at long times (t > T\) During 
the short and long times regimes the film dewetts 
like a Newtonian liquid, with the respective constant ve- 
locities V = |5|/CoA (t < T , A = y/r]oh /Co) and 
Vi = ISI/C0A1 (t > n, Ai = v^iVCo > A ). In the 
intermediate regime (which concerns an important part 
of the experiments, for which the temperature is close 
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to the glass transition temperature) the increase of the 
height of the rim is slowed down by elastic forces, result- 
ing in a linear relation between the width of the rim W 
and the dewetted distance L. The energy balance thus 
reads \S\V ~ (ohoGLV 2 /\S\ which gives a decrease of 
the velocity like . The more rapid decrease (V ~ t^ 1 ) 
observed experimentally Q is obtained by including a 
positive horizontal stress <7o within the film at the be- 
ginning of the dewetting process; the existence of Co is 
attributed to the spin-coating of the PS solution, which is 
followed by a fast evaporation of the solvent, letting the 
PS molecules in a non-equilibrium frozen- in state 

Until now, little attention has been paid to the time 
evolution of the rim width reported in 0, though it is 
most unexpected. Indeed, after logarithmically increas- 
ing, the rim width reaches a maximum, and then slowly 
goes down. One can show that such a decrease of the 
rim width with time is impossible in a classical dewet- 
ting situation where surface tension rounds the rim. In 
the theoretical model described above Q, despite the 
fact that surface tension is ineffective, the rim width of 
a viscoelastic film increases with time during the elas- 
tic regime (W ~ L), while reaching a constant value in 
the long time Newtonian regime: W = Ai = y'jji/io/Co) 
therefore showing no decrease with time of the rim width. 
However, at the sight of the latter expression of W, we 
can anticipate that W could decrease at the condition 
that the friction parameter Co increases with time. 

A non-constant friction parameter is possible in the 
case where the friction force does not increase linearly 
with the sliding velocity. These cases are quite common 
in polymer physics. The friction of a PDMS elastomer 
on grafted or absorbed PDMS surfaces has been shown 
to increase only very slowly with the velocity: Casoli et 
al. found a friction force increasing proportionally to Vs 
on absorbed brushes, and a friction force increasing only 
like Vb on dense grafted brushes, for sliding velocities 
being between lO/im.s^ 1 and 5 mm.s -1 flij . More re- 
cently Bureau et al. showed a V$ dependence of the 
friction force between an elastomer and a grafted brush, 
for sliding velocities ranging from 300/im.s -1 down to 
0.01/im.s -1 which is the velocity range of the dewet- 
ting experiments. 

A general expression of the friction force by surface 
unit is pl| : 

fr = C,V a {£) a (1) 

where a is a shear-thinning exponent smaller than one. 
According to the reported friction experiments, a could 
range between 2/3 and 5/6. The effective friction coeffi- 
cient is then a decreasing function of v: C = Co (V a /v) a , 
where V a is the velocity below which the friction coef- 
ficient remains equal to Co correspondingly to a linear 
friction. Hereafter, we assume V a to be small enough for 
us to omit the linear friction regime. For a Newtonian 
fluid, such a friction force would bring, as an equilibrium 




FIG. 2: Numerical calculations of the evolution of the rim 
width for houo/\S\ — 0, 1, 2, 3, and 4, with a friction exponent 
a = 2/3, and a ratio n/ro = 100. 




100 1000 10000 

Time (s) 



FIG. 3: Width of the rim as a function of dewetting time for 
films (ho = 57nm, M w = 233kg/mol, and T = 125°C) aged 
at room temperature for days indicated in the figure. 

relation between friction and viscous forces within a film 
dewetting toward the x direction: 

, T . ( v \ 1 a d f , dv\ 

^{vj = (2) 

As long as h remains of the order of , the solution of 
Eq.pjl. is simply given by 

v(x,t) = vUl-^-^y (3) 

when x — Vo a t < Ao Q , and v(x, t) = elsewhere. 

Here, the dewetting velocity V 0a = ((^p) V 2 /V a a )^ 
is constant in time, as well as the rim width Ao Q — 

I ((^ir) v o a / y a a )^ A o E3- These analytical results 
have been confirmed by numerically solving Eq. com- 
bined with the volume conservation equation and the 
condition of stress continuity at the front of the film. 

As expected, when the friction force is a non-linear 
function of the velocity, the rim width is a function of the 
amplitude of the dewetting driving forces (i.e. the capil- 
lary forces in usual cases: A 0q ~ |S , |'^ r "). A decrease of 
the rim width can thus be expected if the driving forces 
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decrease during the dewetting process. That is, in fact, 
what happens when a viscoelastic film presents residual 
stresses at the onset of the dewetting. Indeed, the resid- 
ual positive stresses are part of the driving forces, but 
they decrease with time, as the film undertakes internal 
relaxation toward an equilibrated state 0, 0] . 

An equilibrated viscoelastic film would start (t < To) 

dewetting with a velocity Vo a — (Vo/^* 2 ) 2 ""; an d 
a rim width Ao Q ~ (Vo/V a ) 2 ~ a Ao, and then (to < 
t < Ti) slow down to Via — {Vi/V a 2 ) 2= " , ending 
with a rim width Ai a ~ {V\/V a ) 2 ~ a Ai at long times 
(ti < t). When a horizontal positive stress cto is ini- 
tially present, the film starts dewetting with a veloc- 

2 

ity Vi — Voa (1 + hoao/\S\) 2 ~ a and a rim width A, = 
Aoq (1 + h a /\S\)' 2 ~ Tr " . During the elastic regime (to < 
t < ti ) , the residual stress has not yet relaxed much, and 
the rim width increases up to its maximum value: 



h a \ ' 
\S\ J 



(4) 



Around the characteristic time ti , the residual stress has 
decreased significantly, leaving the capillary forces as the 
only driving forces. Consequently, the rim width de- 
creases down to Aia simultaneously with the increase of 
the dewetting distance. Once again, we confirmed these 
results by numerical resolutions of the equations of the 
flow (see Fig. (J3J)). We notice that A m is higher than 
Ai Q only if the exponent a is not nil. The combination of 
non-linear friction together with residual stresses is thus 
necessary to explain the decrease of the rim width with 
time. Interestingly, the residual stresses cro can come off 
from the relation between the initial dewetting velocity 
Vi and the maximum rim width A m : 



Vi = V a 



2(l-a)' 
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(5) 
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FIG. 4: Dewetted distance as a function of dewetting time 
for films {ho = 57nm, M w = 233kg/mol, and T = 125°C) 
aged at room temperature for days indicated in the figure. 
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FIG. 5: Evolution of the maximum rim with as a function of 
the initial dewetting velocity for several aging times and PS 
molecular weights ranging from 35K to 600K. 




The exponent a could thus be deduced from a Log-Log 
plot of Vi as a function of A TO , for dewetting experiments 
with various values of the residual stress. Furthermore, 

we notice a weak dependance of Vi/A„? on the viscos- 
ity 771, and thus on the molecular weight of the polymer, 
when a is close to unity. Then, several Log-Log plots, 
corresponding to different molecular weights, can be su- 
perposed in order to measure a. 

At the same time, we have conducted new exper- 
iments on the aging of thin PS films deposited onto a 
PDMS mono-layer. The consequences of up to three 
hundred days of aging, just below the glass transition, 
on the dewetting have been systematically investigated. 
We observed a slowing down of the dewetting, as well 
as a decrease of the maximum rim width, when increas- 
ing the aging time (see Fig. JHJl and QJ). Eventually, 
at the longest aging time, no maximum rim width was 
observed. We interpreted these results as a manifesta- 
tion of the very slow relaxation of the residual stresses 



FIG. 6: Numerical calculations of the evolution of the dewet- 
ted distance for ho(To/\S\ = 0, 1, 2, 3, and 4, with the friction 
exponent a — 2/3, and a ratio n/ro = 100. 



within the film with aging time. The Log-Log plot of the 
maximum rim with as a function of the initial dewetting 
velocity for several aging times and several PS molecu- 
lar weight revealed a friction exponent a = 0.84 ± 0.1 
(see Fig. (0). The time corresponding to the maximum 
rim width should correspond to the characteristic relax- 
ation time of elastic constraints Ti within the film. We 
notice that this time is significantly shorter than the cor- 
responding relaxation time in bulk {i.e. the bulk repta- 
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FIG. 7: Numerical calculations of the evolution of the dewet- 
ting velocity for ho<To/\S\ = 0, 1, 2, 3, and 4, with the friction 
exponent a = 2/3, and a ratio ti/tq = 100. 



tion time of the polymer chains), which is consistent with 
the fact that the entanglement length should be longer in 
thin films Other structural relaxation mechanisms 

related to the proximity of the glass transition might also 
be involved. 

Another interesting result, which reinforces the impli- 
cation of non- linear friction in dewetting, is that, even- 
though the initial dewetting velocity Vi decreases with 
the aging time, the power-law evolution of the velocity 
during the dewetting process remains the same (V ~ 
i -1 ). This observation is in contradiction with the pre- 
dictions of the theoretical model developed in ||| , where 
a t^ 1 law was predicted for high residual stresses, but 
where the decrease of the velocity should tend toward a 
t~2 law for low residual stresses. On the other hand, 
this observation is coherent with the model we devel- 
oped here. Indeed, the dewetting velocity of an equili- 
brated viscoelastic film goes from Vo a down to V\ a dur- 
ing the elastic regime. The power law of this decrease 
can be obtained from the energy balance, which reads 
\S\V ~ ( h GV a a LV 2 - a /\S\. The later relation gives: 



for To < t < Ti, which is not far from V ~ if a is 
close to unity, as deduced from the plot of W m as a func- 
tion of Vi (see Fig. Q). Thus, the dewetting velocity 
has always a power law close to t , which is made more 
rapid around the relaxation time t\ by the presence of 
residual stresses. This latter result is in complete agree- 
ment with our experiments (see also 0, We have 
also confirmed these analytical prediction by numerical 
solutions (see Fig. (Jfjjl and J7J). 

In conclusion, the non-linearity of the friction force 
between the substrate and the film, together with the 
initial presence of residual stress within the dewetting 
liquid, are necessary to explain the observed decrease of 
the rim width with time. Coincidentally, the non-linear 
friction leads to a rapid decrease of the dewetting velocity 
of a viscoelastic film, which fits the experiments better 
than the velocity decrease predicted by the formers mod- 
els. The good agreement of the predictions on both the 
rim width evolution and the dewetting velocity makes us 
confident about the validity of our model. Finally, the 
Log-Log plot of the initial dewetting velocities as a func- 
tion of the maximum rim width allows to determine the 
exponent a characterizing the friction between the film 
and the substrate. Our study clarifies the role played 
by the friction and the residual stresses, both ingredients 
that may be of great importance in other polymer thin 
film situations. A major challenge for the future would 
be to understand the processes involved in the relaxation 
of the residual stresses below the glass transition. 



V = V U 



(6) 



Acknowledgments 



Pascal Damman is a Research Associate of the Bel- 
gian National Funds for Scientific Research (FNRS). This 
work was supported by the FNRS, the Walloon Region 
and the European Social Fund. 



[1] G. Reiter, J. Forrest, Special Issue on Properties of Thin 
Polymer Films. Eur. Phys. J. E 8, 101 (2002) 

[2] D. G. Bucknall, Progress in Materials Science 49, 713 
(2004) 

[3] F. Saulnier, E. Raphael, and P.-G. de Gennes, Phys. Rev. 
Lett. 88, 196101 (2002); F. Saulnier, E. Raphael, and P- 
G. de Gennes, Phys. Rev. E 66, 061607-(1-12) (2002) 

[4] V. Shenoy and A. Sharma, Phys. Rev. Lett. 88, 236101 
(2002) 

[5] S. Herminghaus, R. Seemann, and K. Jacobs, Phys. Rev. 

Lett. 89, 056101 (2002) 
[6] G. Reiter, Phys. Rev. Lett. 87, 186101 (2001) 
[7] P. Damman, N. Baudelet, and G. Reiter, Phys. Rev. Lett. 

91, 216101 (2003) 



[8] T. Vilmin, and E. Raphael, cond-mat/0502228 (2005). 
[9] F. Brochard-Wyart, G. Debregeas, R. Fondecave, and P. 
Martin, Macromolecules 30, 1211 (1997) 

[10] P.-G. de Gennes, F. Brochard-Wyart and D. Quere, Cap- 
illarity and Wetting Phenomena: Drops, Bubbles, Pearls, 
Waves (Springer, 2003). 

[11] R. B. Bird, R. C. Armstrong, and O. Hassager, Dynamics 
of polymeric liquids (John Wiley & Sons, volume 1, 1977) 

[12] Eventually, at even longer times, surface tension will 
round the rim, leading to a rs decrease of the veloc- 
ity (see 0,0). 

[13] G. Reiter, and P.-G. de Gennes, Eur. Phys. J. E. 6, 25 
(2001) 

[14] A. Casoli, M. Brendle, J. Schultz, P. Auroy, and G. Re- 



■5 



iter, Langmuir 17, 388 (2001) 
[15] L. Bureau, and L. Leger, Langmuir 20, 4523 (2004) 
[16] Even if almost no interpenetration is expected between 
the PS film and the PDMS layer, the non-linear relation 
between friction and velocity could be a consequence of 
the rheologic response of the brush to a shear strain, i. e. 
shear thinning. 

[17] Note that when a tends toward zero, Vja tends toward 



Vo, and Ao Q scales like Ao but diverges. This is a sim- 
ple consequence of the transition from a finite power law 
shape to an exponential shape of the rim. 

[18] L. Si, M. V. Massa, K. Dalnoki-Veress, H. R. Brown, and 
R. A. L. Jones, Phys. Rev. Lett. 94, 127801 (2005) 

[19] G. Reiter, M. Sferrazza, and P. Damman, Eur. Phys. J. 
E 12, 133 (2003) 



